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Using the techniques of out-of-equihbrium field theory, we study the influence on the properties of 
cosmological perturbations generated during inflation on observable scales coming from fluctuations 
corresponding today to scales much bigger than the present Hubble radius. We write the effective 
action for the coarse-grained inflaton perturbations integrating out the sub-horizon modes, which 
manifest themselves as a colored noise and lead to memory effects. Using the simple model of 
a scalar field with cubic self-interactions evolving in a fixed de Sitter background, we evaluate 
the two- and three-point correlation function on observable scales. Our basic procedure shows 
that perturbations do preserve some memory of the super-horizon-scale dynamics, in the form of 
scale-dependent imprints in the statistical moments. In particular, we find a blue tilt of the power- 
spectrum on large scales, in agreement with the recent results of the WMAP collaboration which 
show a suppression of the lower multipoles in the Cosmic Microwave Background anisotropics, and 
a substantial enhancement of the intrinsic non-Gaussianity on large scales. 

PACS numbers: 98.80.Cq,04.62.-|-v 



I. INTRODUCTION 

Stochastic inflation provides an cfflcient approach to 
study inflationary dynamics and has become a very pop- 
ular way to describe the growth of density perturbations 
on scales larger than the Hubble radius. In the first fun- 
damental works P, S 01 S IE 01 1 the inflaton field was 
split into a super-horizon and a sub-horizon part directly 
in the equation of motion. This splitting is operated in 
Fourier space through a window function, that separates 
high from low frequencies. The relevant variable is the 
long- wavelength part, while the sub-horizon modes are 
collected in an effective noise term, playing the role of a 
classical perturbation to the super-horizon dynamics. 

The resulting effective equation of motion is then a 
Langevin-like equation analogous to the one describing 
Brownian motion, where the deterministic evolution is in- 
fluenced and modified by the stochasticity of the source, 
whose effects can be taken into account only as a statis- 
tical average over time. Indeed, in this formalism there 
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is not any knowledge of the exact form of the noise, but 
only of its statistical properties. 

A more general approach exploits the infiuence func- 
tional method H, U, and operates the frequency split- 
ting at the action level getting rid of the high frequen- 
cies via a path-integral over the sub-horizon part of the 
field. The effective action thus obtained contains some 
extra terms that can be interpreted as the coupling of the 
super-horizon field with a classical random noise source, 
whose configurations are statistically weighted by an ap- 
propriate functional probability distribution, becoming 
the origin of the stochastic character of the Langevin- 
like equation of motion. 

The super-horizon degrees of freedom are then treated 
as a purely classical field, all the quantum fiuctuations 
being collected in the classical noise term. This feature 
is claimed not to be a simple computational trick, but 
an intrinsic characteristic of the system. Stochasticity 
is thus not only a clue to understand the properties of 
inflation and the origin of the large-scale structure in 
the Universe, but also as a way to explain the transition 
from a quantum to a classical world From a for- 

mal point of view, the quantum decoherence process in 
the stochastic inflation framework has been discussed in 
various works 0, 0, 0, 0, 0| , where it was pointed 
out that the classicality of the coarse-grained field (im- 
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plicitely assumed in the first papers) is not necessarily 
assured, but is subject to some restrictive conditions. 

Using standard techniques of stochastic processes [T^ 
IT^ , the Langevin equation for the field expectation value 
leads to an evolution equation (the Fokker-Planck equa- 
tion) for its probability distribution function. In the first 
works, the noise correlation time is assumed to be in- 
finitesimally short, and the correlation function for dif- 
ferent times can therefore be considered as being propor- 
tional to Dirac's delta function 6{t — t'), which sets its 
white-noise properties. This assumption allows to apply 
a well-known formalism for the derivation of the Fokker- 
Planck equation and its solution. However, the charac- 
teristic of the correlation function strongly depend on the 
window function, whose choice is not a mere mathemat- 
ical tool, but has several physical effects ^ white 
noise arises only as a consequence of a sharp momentum 
cutoff, whereas a smooth weighting avoids highly singular 
noise correlators and produces a colored noise. 

The choice of a colored noise is interesting for at least 
two reasons. The first is that a sharp momentum-space 
cutoff seems rather unphysical, while a smooth weight- 
ing of the modes is much more likely. Actually, the most 
natural way to integrate out the small-scale fluctuations 
is to average the field in configuration space, choosing 
an appropriate finite volume window function. In most 
cases, this choice results in a smooth weighting in Fourier 
space (thus producing a colored noise), while the sharp 
momentum-space cutoff corresponds to a rather compli- 
cated infinite volume window function in configuration 
space. Moreover, it is possible to single out a wide class 
of window functions for which the shape of the colored 
noise correlation is asymptotically the same A sec- 
ond reason may be the fact that a colored noise could 
play, during inflation, an important role in producing 
intrinsically non-Gaussian density fluctuations as initial 
conditions for the subsequent evolution of the large-scale 
structure of the Universe . 

In the simplest single-field slow-roll models of infla- 
tion, non-Gaussian features in scalar perturbations are 
produced by either the self-interaction of the inflaton 
field j2ll | , which are however constrained to be very small 
by the slow-roll conditions, or by the back-reaction of 
field fluctuations on the background metric, whose am- 
plitude is also strongly constrained by the slow-roll con- 
ditions [23, US 13 123 ■ It been shown, however, that 
the most copious source of large-scale non-Gaussianity 
is stored in the post-inflationary second-order evolution 
of perturbations, which sets in a universal level of non- 
Gaussianity for the gauge-invariant gravitational poten- 
tial, which turns out to be of order unity [26ll27l |. 

In this paper we point out that there is another 
source of intrinsic, and generally scale-dependent non- 
Gaussianity in the fluctuation pattern, which originates 
from the cross-talk between super and sub-horizon scale 
perturbation modes. On scales much larger than the 
Hubble radius, non-Gaussian features generally arise as 
a consequence of the non-linear multiplicative form of 



the Langevi n eq uation, when back-reaction effects are ac- 
counted fo r 1281. 129I I30I [31I ^3 However, as discussed in 
Rcfs. m ilM^Ei ESll El, these effects do not 
directly reflect into the statistical properties of cosmo- 
logical perturbations on sub-horizon scales. Indeed, in 
order to deal with fluctuations relative to our patch of 
the Universe, one cannot simply perform statistical av- 
erages over the entire ensemble of possible states, rather 
one should allow for the observed smoothness of our Uni- 
verse on large scales. A possible, though approximate, 
way to take this constraint into account is to replace en- 
semble averages with averages over the conditional proba- 
bility density functional that fluctuations on sub-horizon 
scales assume a certain value, given that the inflaton field 
is spatially homogeneous at t = fgo, i-e- about 60 e-folds 
before the end of inflation (corresponding to a comov- 
ing scale slightly larger than the present Hubble radius). 
This is equivalent to set, for the probability distribution 
of the fluctuations, the 'initial' condition P(6ip,teQ) = 
5{5ip~Sip6o) ^mM,M,E^M,M,M,E^ Although 
this may happen (and in most models it does) well after 
the beginning of the accelerated expansion, if the noise 
driving the fluctuations is white their evolution is Marko- 
vian, implying that any notion of the previous history is 
erased. The probability distribution then behaves exactly 
as if inflation had started at that time and the level of the 
inflaton non-Gaussianity remains fully negligible. On the 
contrary, a colored noise has a non-vanishing correlation 
time: because of this fact the inflaton keeps memory of 
what happened before the constraint, and its evolution 
ceases to be a Markov process. In this scenario, the prob- 
ability distribution evolves in a different way, and also 
higher moments become important. 

Since the solution of the Fokker-Planck equation with 
colored noise carries several complications and is still a 
partly unknown matter, in this paper we followed a dif- 
ferent approach, trying to perturbativcly determine the 
probability distribution for the inflaton field directly solv- 
ing the Langevin equation in a statistical way. 

The plan of the paper is as follows. In Section |n] 
we briefly describe the derivation of the stochastic equa- 
tion of motion for the inflaton field averaged over super- 
horizon scales, using the influence functional method, 
and evaluate the dependence of the noise on the choice 
of window function. In Section [llll we then choose a spe- 
cific Gaussian shaped filter, obtaining the related (col- 
ored) noise correlation functions, and the variance and 
power-spectrum of the coarse-grained fluctuation field. 
In Section IIVI after introducing a small non-linear (cu- 
bic) term in the potential, we evaluate the bispectrum 
and the third moment of the field. In Section^we inves- 
tigate the memory effects induced by this colored noise 
and build up a formalism to quantitatively determine the 
relevance of the non-Gaussian features of the distribution 
and their sensitivity to the times before tgo. Finally, in 
Section IVII we draw our conclusions. Some technical as- 
pects of our calculation are reported in five Appendices. 
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II. EFFECTIVE SUPER-HORIZON ACTION 

We consider a background de Sitter space-time, whose 
metric reads 



(II.l) 



where a{t) 



is the scale-factor (the Hubble param- 



eter H = d/ais constant in time) and 77 is the conformal 
time defined by drj = dt/a, i.e. rj = —[Ha{t)]^^. 

The action for the inflaton field is the ordinary action 
in curved space-time for a free scalar field with mass m 



(II.2) 



(where greek letters label space-time indices), that with 
our choice of background de Sitter metric becomes 



dxa"^ 
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The equation of motion for such a field is 



3H(j) 



+ m^(f) = 



(II.3) 



(II.4) 



and the standard solution for the inflaton Fourier modes 
reads 
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(II.5) 



where 1^^ — j — jj^ and H'i,^\x) arc Hankcl functions 
of the first kind. In the special case of a masslcss field, 
1/ = 3/2 and 
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(II.6) 



Let us now split the field in two components, divid- 
ing the short-wavclcngth normal modes (with wavelength 
a/k smaller then the horizon scale H^^) from the long- 
wavelength ones. The short-wavelength part is defined 
as 



0±=y dkiy(|k|,t)[(/.k(.T)ak-f (/.^(x)al,]. 



(II.7) 



where the window function W{k, t) projects out the long- 
wavelength modes. The long-wavelength part is then 
simply ^< =(/) — (/)>. Substituting this field decompo- 
sition into the action Ijll.Sp we obtain two distinct free 
actions for the two fields plus an interaction term: 



5[</)<,0>] = 5[0<] + 5[0>] + 5, 



int [V<, 



(II., 



Time Path (CTP) formalism |36l| of out-of-equilibrium 
field theory. Indeed, in ordinary quantum field theory, 
one deals with transition amplitudes in particle reactions 
and one may not study the dynamics of the system. This 
is because one needs the temporal evolution with defi- 
nite initial conditions and not simply the transition am- 
plitudes of particle reactions with fixed initial and final 
conditions. While ordinary quantum field theory yields 
quantum averages of operators evaluated with an m-state 
and an owt-state, the CTP formalism yields quantum 
averages of operators evaluated in the m-state without 
specifying the out-sidXe. In the CTP formalism |33| the 
time integration is made along a closed path going from 
the initial time to positive infinity and back to the ini- 
tial time; the path-integral on the field configurations 
is evaluated on this closed path, along which they need 
not assume the same values on the forward and back- 
ward branches of the time contour. This is equivalent 
to considering two fields <j)^ and (/)^ , with the constraint 
(/)^(+oo) = (/)^(+oo), with the ordinary single time in- 
tegration on the real axis. These two fields have to be 
set equal to each other in the equation of motion. We thus 
end up with four different fields 0+ , (/) J , (/)~ , (/)~ for small 
and large scales, and forward and backward times. 

The effective equation of motion for the in-in expecta- 
tion value of the super-horizon fields 0^ = ip^ can be de- 
rived by integrating the action over the sub-horizon field 
variable 0> . The super- horizon effective action obtained 
in this way contains two extra terms in addition to the 
ordinary action ljll.3|) . describing physical effects related 
to the horizon crossing of the various normal modes 0. 
One of them contains both dissipation and non-local mass 
renormalization effects, while the other describes the in- 
ffuence on the super-horizon modes by the sub-horizon 
ones, whose quantum fluctuation can be treated as a 
stochastic noise. This second term is purely imaginary, 
and as such it cannot be interpreted as a standard action: 
indeed, it appears as the result of a statistical weighting 
over the configurations of the stochastic noise fields repre- 
senting sub-horizon quantum fluctuations, which couple 
to in the effective action S'off[<p^]. 

After some manipulations (see Appendix ^ for de- 
tails), introducing two real classical fields ^1 c ^2, whose 
configurations are statistically weighted by the probabil- 
ity distribution functional 7^[^i(a::), ^2(2;)], it is actually 
possible to write 



(e'^"''[^*l)s = Jv^iV^2n^ui2y'^''^^^\ (II.9) 



In order to obtain real in-in vacuum expectation values 
for the field (j) instead of the usual in- out transition am- 
plitudes, the standard procedure is to work in the Closed 



where the effective action S'off[<p^] is (introducing ip^ 

tp+ - if-) 
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S^s = S[ip+] - Slip"] + / d'^xaiH 



H 



(11.10) 



The quantum noise on small scales acts then as a classical random source. Here, S is simply the free action (|II.3|I . 
since the extra dissipation and mass renormalization terms are small and can be neglected on a first approach, as 
shown in Appendix IbI 

The statistical weight of the two random fields and ^2 is Gaussian, 



(11.11) 



where A ^(x,x') is the functional inverse of the symmetric (under simultaneous exchange of discrete and continuous 
indices i,t j, t') kernel 



A(x,x') = — /— — k^{7]7]'f^^W'{kj])W'ikr]')Rc[M,{kr],kT]')], 

STT I k kr 



(11.12) 



with 



M^{kri,kT]') = 



Hi'\k\r^\)Hi'>{k\i\) -kiHi}\kM)Hl'2l{k\r^'\) 
^kr^Hl'},ik\rj\)Hl'>ikW\) k^\rjMHl'2,{k\7^\)Hi'}likW\) 



(11.13) 



In the procedure we described, quantum fluctuations of 
the sub-horizon modes are collected via the path-integral 
in a rapidly varying classical noise term coupled to the 
super-horizon part of the scalar field: these fluctuations 
can thus talk to the super-horizon modes and perturb the 
dynamics on scales larger than the Hubble radius during 
inflation. 

The coupling of the two random noises ^1 and ^2 with 
the scalar field is slightly different from that in Ref. 0; 
the reason is that the choice of a general window function 
W{kr]) (which is not necessarily able to produce a sharp 
cut in the frequencies, but can have some spread around 
the horizon scale) can introduce a fc-dependence in the 
effective field coupling to the noise, thereby spoiling the 
separation between super- and sub-horizon scales. Our 
different choice of basis avoids this problem. Another 
consequence of this formulation is that the correlation 
functions of the two noise fields are different from each 
other , an d cross-correlations appear, in a similar way to 
Ref. [11. 

The effective equation of motion obtained from this 
action with the usual CTP method is then 




= (/3 + 3H(p 



-^+m\^3H^i-ii+H^2, (H.14) 

that is a Langevin-like stochastic equation, where the 
dynamics of the field (p is subjected to random "kicks" 
given by the rapidly varying stochastic force ^. We treat 
the effect of the random force as a perturbation of the 



classical dynamics and split the field if into its mean, 
obeying the classical equation of motion (jll.4p . plus a 
fluctuation Sip that by definition vanishes when averaged 
over all noise configurations. In the massless case, the 
equation for the fiuctuations becomes 

Sip + 3H6ip - = 3i7a + 6 - H^2] (H.15) 

if we neglect the exponentially suppressed spatial gradi- 
ents, the second time derivative dip (assuming the validity 
of the slow-roll conditions) and the ^1 term, we finally get 

S^ = ^i^j^C (11.16) 

Even in this simple form, a deterministic treatment of 
this equation is impossible, since we do not know the 
exact form of ^1 and ^2- However, we can study this 
equation from a stochastic point of view [l^ IT^ , in or- 
der to understand how the statistical properties of the 
Gaussian noise (that are completely characterized by the 
two-point correlation functions Aij{x,x')) determine the 
behaviour of Sip[£], now treated as a stochastic variable 
itself. In other words, our goal will not be the exact de- 
termination of the evolution of the field dip^ but of its 
probability distribution functional. 

III. SPECTRA AND TWO-POINT 
CORRELATION FUNCTIONS 



From the equation of motion Ijll.lGp . one can immedi- 
ately compute the two-point correlation function for Sip, 
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which reads 

{S^{x)6^{x')) = / di/'k'(^(t,x)e(F,x')) , (III.l) 

teo'^ tea 



where the correlation function of the noise ^ is given by 



1 

1p 



Aii{x,x') 



Ai2{x,x') 



A2l{x,x') A22{X,X') 
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(III.2) 



and the exact form of the matrix elements Aij depends 
through (|11.12() on the choice of the window function 
W{kr]). 

It is interesting to note that if we project the modes 
using the Heaviside step-function W{h]) = d{{k/aH) — 
e), in such a way that W'{kr]) = S((k/aH)—e) = aHS{k— 
eaH), the noise correlator gives the standard result 

obtained in the first stochastic inflation works 111, 



sin eatHr 



{l + e^)S{t-t'), (III.3) 



while the other correlators are of order 0{e^). Therefore, 
in the small e limit, we get the standard result (for x = x' 
and t < t') 



{S^{t)5ipit'))^^{t-t,o). (III.4) 



Integrating out long wavelenghts using the step- 
function gives what in stochastic language is called a 
Markov process. However, this is not the most natural 
choice one can do. Actually, the smoothing of the field 
is generally performed in configuration space through a 
fimction w{r/R) that rapidly decays for distances much 
larger than R. In momentum space, this operation pro- 
duces a weighting of the modes with the Fourier trans- 
form w{kR) that projects out the high frequencies. Since 
we are interested in the short wavelenghts part of the 
field, and our smoothing scale is the comoving Hubble 
length 1 77 1, a natural choice of the momentum window 
fimction will then be W{kr]) = 1 — w{k\rj\). 

In a general case the choice of a more physical way to 
separate the modes gives a colored noise term, which is 
unfortunately much more difficult to treat. Namely, if we 
smooth the field with a Gaussian filter 



'w{r/\ri\) = e 



(nL5) 



we get the window function 



W{kr^) 



(in.6) 



also in this case we can compute the noise correlation 
function for v = 3/2, and in the limit r = |x — x'| ^ 
we obtain (Appendix [nj, setting t = t ^ t' , 



where 



47r2 cosh^ Ht f-^""'^ 

A.;— 



>^^(fe)(^)(-l)'=(fc + l)(fc + 2) 



(2fc- 1)!! 



sinh^ 



cosh Ht 



(111.7) 
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cosh Ht 



cosh Ht 
2k 
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2k 



cosh Ht \ ' e^^ 
4 (fc + 3) 



- 1 



cosh Ht 
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For generic values of the parameter cr, the noise corre- 
lation function (^(i)^(i')) '^^^ have a rather complicated 
functional form, and is plotted in Figure lllll However, in 



the small a limit, the leading term of the scries (jIII.7|) is 
that with k = 0. Therefore, all the Aij vanish but An, 
and the noise correlation function becomes 1X911 



(mm) 



i/4 



1 



87r2 cosh^ {H{t-t')) 



0{a^), (IH.9) 



tion gives (to leading order in ct^) 



cosh(g(t - ^6o)) cosh(g(^^ ~ teo)) 
^ cosh{H{t -t')) ^ { ■ ) 



which inserted in (jllLip and after a double time Integra- that for t' ^ t ^ tgo approaches the standard result 
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Figure 1: Noise correlation functions {(.{t)^{t')) for r = and 
five values of a varying from 1 (top line) to (bottom line). 
The only non- vanishing function for cr = is (CiCi) 



Conversely, for i = t' tgo one obtains 



— r In cosh" {H{t 

OTT 



^4 
OTT 

(III.ll) 

we thus find that soon after the initial condition is set 
the variance depends quadratically on time, unlike the 
white-noise case (|III.4|I . where the time dependence is 
linear. This behaviour, which comes from the fact that 
the Dirac delta function is divergent while a more phys- 
ical correlation function is not, will have an important 
role in the following analysis. 

The power-spectrum can be immediately derived by 
inserting the explicit form l|II.12() of the Ai 



into (HTLTJ. 

The two integrations over time factor out, and per- 
forming the change of variables t x = —kfj and 
t' ^ y = —kfj' we get 



For 1/ = 3/2 and using the Gaussian window 1III.6|I . so 
that W'{x) = a;exp[— fc^77^/2a-^]/cr^, the functions /i and 
/2 are given by 



/i(x) = 5 e , 



f2{x) 



3a2 



re 2o 



(in. 16) 
(in.17) 



Setting t = f' va (|111.12|l . we extract the power-spectrum: 



2ti J 



k\ri6o\ , 

dx(fi{x) - f2{x) 

k\v\ ^ 



(in.18) 



For t 3> ^60 J we can take into account scales that are 
much larger than the horizon at time t (such that —krj <Si 
1) but much smaller at tgo (— fc'/60 S> 1); for such scales 
we get 



dxfi{x) ~ dxf2{x) 



(in.19) 



Studying the behaviour of the integrals of /i and /2 for 
large and small values of the parameter a we see that 
for cr ^ wc have /g dxfi{x) i and /g dxf2{x) — > 0, 
while for a ^ oo both integrals tend to vanish. This 
means that for small a the power-spectrum on very large 
scales (such that k <ti aaH) approximates the standard 
result i/^/47r^, while fluctuations do not appear at all for 
very large values of the parameter. This behaviour is due 
to the fact that for small a the window function tends 
to unity, and larger and larger scales are included in the 
noise ^, whose correlation function therefore reproduces 
the ordinary fluctuation behaviour (when all scales are 
taken into account). For increasing a, instead, since the 
window function vanishes, the noise contains only smaller 
and smaller wavelengths, and is no more able to influence 
super-horizon scales. 



{5ip{x)S^ix')) = 



dk sin kr 
k kr 47r^ 



F{kri,kT]'), (IIL12) 



with 



F{kT], kr]') = Re 



and 



/i(^) 



k\rieo\ 

dx{h{x) - h{x)) 
X / dy{h{y) ~ f2{y)y 

Jk\v'\ 



^x'/'W\x)Hl'Hx), 



, (in.13) 



.5/2 

2 "3~ 



W'{x)hI'Ux). 



(in.14) 
(in.15) 



IV. THE INTERACTING SCALAR FIELD 

Let us now introduce a small self-interaction term in 
the potential for if, in the form of a cubic term 



V 



(IV.l) 



where 



<^ 1. This toy- model may be useful even to 
describe the generation of non-linearities in scalar fields 
other than the inflaton field and then trasmitted to the 
latter, as described in Refs. [sM l39| . We expand to 
second order in the fluctuations the field Sep substituting 
Sip = Sifii + Sip2 in Ijll.lBp . Assuming that the small self- 
interaction term can only produce second-order effects 
(acting therefore only on 5ip2), we can solve recursively 
the equation of motion for Sipi and 6ip2- In the slow-roll 
approximation they become 



3H6ipi = 3H£_i + a 
3HSp2 + ^Sifl = 0. 



(IV.2) 
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To first order in the field expansion, the three-point func- 
tion decomposes into the sum of three terms: 

+ {Sip2ixi)S(pi{x2)S(pi{x3)). (IV.3) 

Since Sip2 depends quadraticaUy on Sipi, each of these 
first-order terms is actuaUy constituted by the sum of 
three four-point functions of 6(pi, two connected and one 
disconnected. However, since in the perturbative expan- 
sion in of the equation of motion the originally vanish- 
ing mean value of Sip is shifted away from zero by {Sip2), 
and the disconnected term is indeed proportional to this 
quantity, we can eliminate it just by requiring that the 
mean of the fluctuation vanishes. For this purpose, we 
set 

= f dt'[5vl[t',^) - (<5^?(t',x))]; (IV.4) 

the second term, added in order to make the mean value 
vanish, cancels the disconnected contribution, while the 
connected ones are equal to each other. The first term of 
(|IV.3p . namely, becomes 



where 



dt'(<5^i(ti,xi)(5^i(t',X3)) 

X (,5(pi(i2,X2)<5(pi(t',X3)), (IV.5) 

that is (setting ii = t2 = ^3 = i and going to conformal 
time) 



3i?2 /^^^^ ,/ 4(27r)6 J 
fdk' 



X f^e'^'(^^-^^)F{k'ri,k'7]'). (IV.6) 

Going to Fourier space and extracting a factor of 
(27r)^(5(ki + k2 + ka), we obtain the total bispectrum 
as the sum of the three contributions 



-B(ki,k2,k3) = Bi2 + permutations. 



(IV.7) 



gix) = 1 + y 



1 + 
6 



(IV.IO) 



and the contributions to the total bispectrum can be 
now calculated analytically. In the limit of super-horizon 
scales {—kiTi/cr <^ 1, but — fcijyeo/c S> 1) and to order 
0((T^) we obtain 



'^l'^2 
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2cr2 
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(IV.ll) 



If we take as our window the step-function, so that the 
derivative W'{x) is Dirac's delta function (5(.t — e), we get 



F(fc77, fc,7') = ( 1 + y + y ) 



fcl^l) 

79(e-fc|77'|)z^(fc|7,6o|-e) 



and 

Bl2 



fiH^ 1 

12 A; -j^ 



1 
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(IV.12) 



(IV.13) 



with k^axM < e < fc,„.„,l??ml. This result is in good 
agreement with the one in Ref. |2l| . the main difference 
being that the relevant scale is not ki + k2 + but kmax , 
which is due to our use of the slow-roll approximation, 
instead of the exact solution of the equation of motion. 

For later convenience, let us also calculate here the 
skewness of our coarse grained inflaton field: from ljIV.3|) 
and ljIV.5p . with Xi = X2 = x^, we have in the white- 
noise case (when {S'pP') oc i — igo) 



(V(t,x)) 



(IV. 14) 



while with our window function the result, in the small 
(T limit, is 



(V(t,x)) 



4(27r)^ 



In 



iv^ + VmW + Tim) 



(IV. 15) 



where 



Bl2 — 



" dr]' 



12 /cjfcl J,^gg rj' 

X F{ki7^Mv)F{k2Tl,k2ri'). (IV.8) 

With good approximation, we can write for the Gaussian 
window case 



for early times it can be approximated by 



(V(i,x)) 



12 (27r)' 



(IV.16) 



while for t ^ igo we asymptotically recover the ordinary 
resuh IjlV.Up . 



F(krj, kr]') 



e 2a 



~gikr]) 



X e gikrf) — e 



'g{kmo) 



.g(fc776o) +0(a4) (IV.9) 



V. CONDITIONAL PROBABILITY 
DISTRIBUTION AND MOMENTS 

In this section we derive an expression for the proba- 
bility distribution of the field 5lp at time t as a function of 
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the stochastic variable C = Ci ~ ^2/3. More precisely, we 
are interested in the conditional probability Pt{Sip\SipQQ) 
that the stochastic variable S(p[^]{t) assumes a specific 
value Sip at time t given that it assumed the value Sipea 
at the earlier time <60i '-e. 60 e- folds before inflation 
ends. In particular, we will consider dpeo = 0, as ar- 
gued in 1^ 13 This is the crucial quantity that 
distinguishes Markovian processes from non-Markovian 
ones. A stochastic process is said to be Markovian if its 
conditional probability depends only on the value of the 
variable at the time when we put the constraint, while 
in non-Markovian cases a memory of what happened at 
earlier times is kept. 

It is a well known result of the theory of stochasticity 
jlfii Il7l | that Langevin equations with white noise de- 
scribe Markovian processes, while colored noise is a typ- 
ical cause for Markovianity to break down. Thus, in our 
case the fluctuations of the inflaton field on large scale 
behave like a Markovian process if we choose W{kr]) to 



be a step-function, that generates white noise, but be- 
come non-Markovian whenever a smoother separator is 
applied, as for example our Gaussian window function. 

Our goal is to estimate the amount of non-Gaussianity 
produced in this latter case: in the standard scenario, 
non-Gaussian features are generally small, since they do 
not have enough time to develop after tgo , when the con- 
dition SpQQ = applies. However, if the constraint does 
not erase the memory of earlier times, non-Gaussianity 
could be significantly larger. 

According to Bayes theorem, the conditional proba- 
bility is obtained from the joint probability Pt {Sip, Spqq) 
that the random process S(p[^] assumes the values Speo 
at time igo and dp at time t, normalised with the prob- 
ability P{6ipeo) that Stp[^]{teo) = Sip^Q. A simple way to 
evaluate the joint probability (see e.g. Rcf. |0) is to 
average over all ^ configurations the product of two delta 
functions centered on these values: 



Pt{S^,6p^eQ) ^Mlvi6(6p>m)~5p) 5(5pmUo) - 5p>Qo 

I 



(V.l) 



where the shorthand notation ^-^A^-"^^ stands for the 
double integration in the exponent of the Gaussian 
The probability distribution P{5lpqq) 



weight ljll.ll|) . The probability distribution P((5(p6o) for 
5lp[^\ at teo can be obtained in the same way by averaging 
just one delta function: 

(V.2) 

The equation of motion for 5p[£] is a first order dif- 
ferential equation, and therefore for a given choice of the 
function ^(t) the solution depends on the initial condi- 
tion 5pin at i = tin defined to be time at beginning of 
infiation. 

Since 5p) is a fluctuation (i.e. with vanishing mean 
value), we require that the solution vanishes at every time 
when averaged over all ^'s: the initial condition must 
then be Sipin = for continuity. 

If we could solve the equation of motion, we could then 
calculate the value of the fluctuation at <60- However, 
because of the presence of the delta function, this value is 
constrained and the path-integral runs only over the noise 
configurations that satisfy the condition on 5p>[^]{tQQ). 
Moreover, the value 6p(,Q assumed at tgo constitutes a 
new initial condition for the equation of motion at later 



times, whose solution thus depends only on the noise 
configuration after tgo- 

Writing the equation of motion in integral form and 
then solving perturbatively by iteration, to first order in 

-m ""^ 8*5* 



*60 ,, ptQQ 



«(t6o) = y^ dt'm--^j^ dt'sp'mt') (v.a) 

that depends on the noise term configuration only up to 
^60 5 a-nd for later times 



Spm)^Speo+ / dt'at') 



6H 



Jteo 



JL 

6H J, 



dt'f fdt"at")) , (V.4) 

y'^ieo / 



where in the last equality we did not include the first 
order correction to the linear term in ^ because its effect 
would just be a sub-leading contribution to the two-point 
correlation function. 

As we already pointed out, the solution after igo in- 
volves only integrals over later times: the only possibility 
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for the fluctuation field to keep memory of earlier times is 
then that this configuration itself is influenced by the con- 
figurations before the constraint, which is exactly what 
happens in the case of colored noise. 

We see that in these equations there are both linear 
and quadratic terms in ^: the quadratic terms can be 
added to ^^A~^^ of Eq. IjV.f p in such a way as to obtain 
a modified integration kernel and perform the Gaussian 
integration over the ^'s. Skipping all the technicalities 
(see AppendixEJ), after inverting perturbatively the new 
integration kernel we obtain for the conditional probabil- 
ity the following result: 



Pt{6ip\Sifi6o) = 



2n /do;, e~*'^'^="^e"^'^^^"2+^67T'D222ai 



(V.5) 



where we adopted a convenient short-hand tensor nota- 
tion, setting 



5ipi =6ip~ SipQo + -^^^loit " *6o) 



S(p2 = Sip 



(V.6) 



60 



c„ = jdt'jdt"{amt")) , (V.7) 
A,.= idi fdt'U'iamn) ldT'(dT"{ar')ar")) 



li Ji 



Ik J I 



(V.8) 



and Ii — [<60,^]i I2 — [timtw]: h — [^60i^- I2 — [timt] 

for the integration supports. 



For white-noise processes, the only terms that survive 
are Cn, C22, Dm and -D2221 while those with mixed 
indices vanish. For instance, we have 



C12 (X I dt' dt"5{t' - t") = 0, (V.9) 

J tea Jti-n 



tea / rt rt ^ 2 



Diaioc / dt{ / dt' dt"6{t' -t")\ =0 (V.IO) 

■Jtir, \Jteo Jti„ J 

since the integration supports of the two time variables 
are disjoint, and the same happens for all other mixed 
terms. In this case, thus, the two integrations over ai 
and a2 in the numerator factor out, and the second one 
simplifies with the one in the denominator. We then get 



Pt{6ip\Sip6o' 



2ttJ 



exp 



6i? -'^11195^ 



exp 



2Cii 



(V.ll) 



i.e. a probability distribution function for a Markovian 
process, that does not keep any memory of what hap- 
pened before the constraint (indeed, this function con- 
tains only integrations over times after tgOj since no in- 
dex 2 appears). To zeroth order in this distribution 
is Gaussian with mean Sipeo and variance Cn, which is 
exactly the standard case of equation Ijlll.ip . 

In a case with colored noise, the correlation function 
is no more a delta function but has a finite width, and 
the coefficients with mixed indices no longer vanish even 
if the two time integration supports are disjoint. In the 
general non-Markovian case no factorization is possible, 
and after integrating over ai and a2 we can write 



Pt{Sip\5ipea) 



1 exp 


\ D--1 1 


exp[-i[C %Sip,Sipj] 


y/2'KCn{l-y) exp 


6ff 222 35^3 


exp 


2 C22 





(V.12) 



where C is the matrix with elements C,, and 



1 



y 



C12 
1 

C22 



y ■ 



^2 



C11C22 



(V.13) 



It is immediate to check that from this general formula we can get as a particular case the Markovian formula by 
simply setting C12 = [i.e. y = 0) and Ajfe = for {ijk} ^ {III}, {222}. 

Expanding the derivation operators to first order in ^ , after some algebra we get 



exp 



Pt{6ip\Sipeo) 



2 Cii(l-y) 



\/27rCii(l-y) 
/ 



V 



6H 



D 



D 



222 



dSip 



2C22 



e 2C22 



(V.14) 
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To zeroth order in ^ we then have again a Gaussian 



6H 

distribution, but now the mean is 



C 



{S(p)\nm = Sipea { 1 + ^ ) , 



(V.15) 



that is the sum of the constraint value plus an extra term 
due to non-Markovian memory effects, while the variance 
becomes 



{S^^)\nm = Cn{l-y). 



(V.16) 



In order to calculate the skewness, we have to take into 
account all the first order contributions. The full calcu- 
lation gives 



{Sip^)\NM 



+00 



dSif Sip^ Pt{6ip\SipGo) 



-^Dni{l + K), (V.17) 



where K is given by 



K 



C12 2Dii2 + D121 Cii 1 - 3y 2Dr22 + ^212 



22 
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^22 2 Dill 

C11C12 3 - 5y D222 



{C22f 



D 



111 



(V.18) 



y and K are then the terms monitoring the importance of 
memory effects induced by the colored noise in the con- 
ditional probability distribution variance and skewness 
respectively. Whenever these coefficients are small (as 
for white-noise cases, when they identically vanish since 
C12 = and D112 = D121 = -D122 = D212 = 0) we get fo r 
the variance {5lp^)\m = Cn the ordinary result (|III.1|1 . 
while Eq. ljV.17|) reduces to 



dt 

^60 



M'fdt"{amt"))] 

30 ^iQO 



(V.19) 



which is exactly what we obtain from Eq. (|IV.3p - ljIV.5|) 
with xi = X2 = X3. If instead y and K are significantly 
different from zero the effect is big, and the procedure 
of neglecting times before tgOi a-s we did to derive for 
example Eqs. Ijlll.lip and (|IV.15p . is not correct any 
more. 

We now want to apply the formalism developed so far 
to the case of the Gaussian window (|III.6|I . It is a generic 
feature of non-Markovian systems that, since memory ef- 
fects appear through secular terms which depend on the 
whole time interval before the constraint, they increase 
as this interval increases. Conversely, much after the con- 
straint these effects tend to be erased as a consequence 
of stochasticity. The results derived in Section UTTl should 
thus still hold for t — tgo ^ H^^, while we expect the 
difference from the Markovian case to be maximal for 
t-Uo-^ H-^ and ieo - Un > . 

In this limit, we can compare the white- noise corre- 
lation functions (|III.4|) and (|IV.14p with the ones just 



derived. For the two-point function, while in the Marko- 
vian case we had {Sip^) ex H(t — tQo), now, since the noise 
correlation never diverges and Cn contains a double inte- 
gration over time, for t tgo we get {S(p'^) oc tgo)^, 
and the ratio of the two vanishes. Precisely, for a value 
of the parameter a not too close to 1, at the Cij read 



Cii 



27r/ 2r]T]Qo 



47r2 2 



(V.20) 



Ci 



C'21 



2 (7?2 



47r2 2 



(V.21) 



C22 — 



In 



47J-2 



[H{teo-U„)^ln2]. (V.22) 



We thus get for the y parameter 



C11C22 2H{tQQ~tin)^ 



(V.23) 



and the variance becomes 



1 



Sn^y 2H{teQ-U„), 

(V.24) 

This expression shows two distinct effects. The first is 
the fact that the variance is an increasing function of 
the total amount of e-foldings of inflation before the con- 
straint, and since y tends to vanish (though not very fast) 
for teo — tin ^ its dependence saturates to the value 
of the unconditional variance (jllLlip . This is a physical 
consequence of the fact that, even though highly oscillat- 
ing noise configurations arc statistically suppressed due 
to the non-vanishing correlation time, if there is much 
time after the beginning of inflation the noise about teo 
is almost uncorrelated with its initial value, while if the 
starting time is closer to teo there cannot be very large 
fluctuations of the noise ^, and also the variance of Sip 
gets smaller. A second effect is the quadratic time de- 
pendence soon after the constraint, versus the linear de- 
pendence of the standard Markovian case (|III.4p . This 
behaviour, already seen in Eq. 



(|III.11|I . does not have 
anything to do with memory effects, but is merely due to 
the local shape of the noise correlation function, which is 
never divergent. Actually, we still have that 



(V(t))|jVM AiV 
(VW>Im ^ 2 ' 



(V.25) 



where we have defined AiV H{t-teo) = 60-iV, with N 
the number of e-folds till the end of inflation (see the inset 
in Figure |2J|. This result shows that the power-spectrum 
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Figure 2: Left panel: plot of the variance {Sip'^{t)) as a function of AA'^ — H{t — teo), in the non-Markovian (thin continuous 
lines) and Markovian (thick dashed line) case. Right panel: plot vs. AA*' of the ratio ^^1^^,2(1)1"" between the non-Markovian 
(NM) and Markovian (Af) case. The inset contains a magnification of the region enclosed in the dashed box, showing in 
detail the behaviour for small AA^ (i.e. close to teo)- In all plots, different continuous curves represent values of H{t6o — tin) 
corresponding to 1.5, 2.5 and 5 (from bottom to top). The conditional variance is at any time after teo an increasing function 
of the number of e- folds before tao, converging to the value of the unconditional variance 1111. lH 




0.01 0.02 0.03 0.04 0.05 12 3 4 5 
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Figure 3: Top panels: plot of the K coefficient for the Gaussian window vs. AA^, showing memory effects in the third moment 
{Sip'^{t)), over two different time scales. Bottom panels: dependence on AA'^ of the third moment of Sip (left panel), in the 

non-Markovian (thin continuous lines) and Markovian (thick dashed line) case; behaviour of the ratio (right panel) 

between the non-Markovian (NAI) and Markovian (M) case. Different curves correspond to values of //(teo — ti„) from 5 
(bottom line) to 30 (top line). 
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of perturbations on large scales is naturally bluer than 
the standard one. This basic conclusion is in qualitative 
agreement with the recent WMAP results which show a 
suppression of the lower multipoles in the Cosmic Mi- 
crowave Background (CMB) anisotropics. 

While memory effects do not clearly show up in the 
variance, for the third moment we obtain, since K is dom- 
inated by the D222 term (see the explicit computation of 
all the Dijk coefficients for this choice of the window in 
Appendix IE)| . 



3/,\\i _ A* 3C11C12 
{dip [t))\NM ^ -— - — D-}. 



H 2C22 



^in ) 



;(2^r 



and therefore (recalling (|1V.14|| ') 

(VW)Ua/ 3 



H{t(if) — ti 



(V.26) 



(V.27) 



{S^HmM 8 

showing that for wavelengths slightly smaller than the 
present horizon, the third moment can be considerably 
larger than it is usually assumed, depending on the pre- 
vious duration of inflation. However, at later times (for 
AA^ ^ 1) both the variance and the third moment con- 
verge to the standard Markovian behaviour. The precise 
value of all these quantities for generic times is shown in 
Figure 121 and Figure O 

The amount of non-Gaussianity can be evaluated from 
the ratio of the three-point function and the variance to 
the appropriate power. It is useful to consider both the 

quantities R3/2 = (gy?(t))3/2 and R2 = (s^'ilt))i (the sec- 
ond one is proportional to the non-Gaussianity strength 
parameter Jnl, see e.g. Refs. I^HS) as a function of 
H(t — teo) and Hit^Q — tm). In this case too, as we see in 
Figure^ the effect increases as the time t gets closer to 
teo and farer from tin ■ Analytically, in this limit we have 

M H{tQQ — tin) 



R 



(NM) 
'3/2 



{NM) 



H 4x/2^ ' 

A* H{tQO — tin) 



i?2 



2AN 



(V.28) 
(V.29) 



(where the superscript NM stands for non-Markovian) 
while for the white-noise case 



r: 



(M) _ M (A^)'^' 



3/2 



H 6tt 



R. 



{M) 



AA^ 



3 



(V.30) 
(V.31) 



(where the superscript M stands for Markovian). In 
the Markovian case both quantities vanish, but in non- 
Markovian ones they are finite or even divergent for 
AA^ — !■ 0, because of the quadratic (instead of finear) 
time dependence of the variance, and are growing func- 
tions of H{teo — tin), due both to the increase of the third 
moment and the decrease of the variance with the total 
number of e- folds before t^Q. Non-Gaussianity can then 
be larger by orders of magnitude in our scheme. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we used the stochastic inflation approach 
to address an important theoretical issue, namely that 
of estimating the possible influence of super-horizon per- 
turbation modes on the statistics of cosmological pertur- 
bations on observable sub- horizon scales. To this aim we 
flrst modifled the standard scheme in order to allow for 
the cross-talk of perturbations on super and sub-horizon 
scales. Such an effect is indeed totally absent in the 
traditional stochastic inflation dynamics; this is an ar- 
tifact of the sharp k-space filter adopted in the coarse- 
graining procedure, which is easily removed by adopting 
a smoother filter function, leading to a colored - rather 
than white - noise source in the Langevin-like equation 
which governs the evolution of the coarse-grained inflaton 
field. This modification implies that the evolution of the 
coarse-grained inflaton field behaves as a non-Markovian 
stochastic process, in contrast to the standard case. Per- 
turbations relevant to our smooth local patch of the Uni- 
verse are then consistently defined by constraining the in- 
flaton field to be homogeneous at a conventional time fgo 
(i.e. about 60 e-folds before infiation ends), correspond- 
ing to the horizon crossing of a scale slightly larger than 
the present Hubble radius. As a consequence of its non- 
Markovianity, the coarse-grained inflaton field preserves 
some memory of its dynamics prior to the constraint, 
which means that sub-horizon-scale perturbations have 
some knowledge of the state of the Universe on super- 
horizon scales. 

Endowed with this extended stochastic inflation scheme 
we are finally able to calculate the conditional second- 
and third-order moments of inflaton perturbations on ob- 
servable scales. We perform our calculations in the case 
of a simple model where a scalar field with a small cubic 
self-interaction term evolves in a fixed de Sitter back- 
ground. Our most important findings are: 

• The variance of inflaton fluctuations grows 
quadratically with time around igo- and is there- 
fore smaller than in the standard Markovian case, 
which is linearly dependent on time; this is equiv- 
alent to say that the power-spectrum of den- 
sity perturbations gets bluer on very large scales, 
without invoking any ad-hoc new physical input. 
This generic feature looks very intriguing, since 
the CMB anisotropy power on the largest angular 
scales observed by WMAP appears to be lower 
than the one predicted by the standard model of 
cosmology with almost scale-invariant primordial 
perturbations arising from a period of infiation. 

• The skewness of inflaton fluctuations is larger than 
the standard case at times around teo, which 
is equivalent to an enhancement of the non- 
Gaussianity level on large scales, for a given value 
of the self-coupling strength. The possible pres- 
ence of large-scale non-Gaussianity is an important 
prediction to be confronted with current observa- 
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Figure 4: Top: Behaviour of -R3/2 (left panel) and R2 (right panel) as a function of AA*', for the non-Markovian (thin continuous 
lines) and Markovian (thick dashed line) case. Bottom: behaviour of the ratios -^3/2*^ (^^^ panel) and b!'^'^'^ / R':2 '^ 
(right panel) of the two cases. In each plot, different continuous curves represent values of H{teo — tin) varying from 5 (bottom 
line) to 30 (top line). 



tional limits |42|. Let us also mention that some 
recent analyses have reported evidence for a pos- 
itive detection of non-Gaussianity in the WMAP 
1-year data [H H El E 113 ■ 

In spite of the simplicity of the considered model, we 
can conclude that the cross-talk between super- and sub- 
horizon-scale perturbations is an important effect which 
would deserve an accurate treatment. The scheme pre- 
sented here should be considered only as a first step in 
this direction, which needs to be largely improved and 
extended in various ways. First of all, one should consis- 
tently include metric perturbations in the inflaton dy- 
namics, going beyond our fixed de Sitter background 
treatment. Second, our choice of a Gaussian filter can 
be somewhat arbitrary, and it would be important to 
understand how much the results depend on its choice. 
The conclusions reached in Ref . , according to which 
there is a wide class of filters whose noise correlators show 
the same asymptotic behaviour {S,{t)S,{t')) oc e~^-^l*~* I 
found in this paper (see Eq. Ijin.9|l ). are however quite 
encouraging in this respect. Actually, since the window 



function is only a technical device, we would expect our 
results to be independent of its choice, at least qualita- 
tively. The white-noise exception is not significant, being 
merely a consequence of the "bad" choice of smoothing 
in configuration space. 

The existence of memory effects means that the de- 
tailed dynamics of inflation plays a role in the specific 
form assumed by observable quantities. For instance, in 
our model one finds a residual dependence on the time 
when inflation started, ti„. Nonetheless, it is quite likely 
that this dependence saturates to a universal value deter- 
mined by the general asymptotic behaviour, if the overall 
number of inflation e-folds is very large, as it is usually 
the case. More difficult is to avoid a dependence on the 
precise time tgo sti which we put the homogeneity con- 
straint, and on the specific form of the constraint, which 
we introduced here through a delta function in field con- 
figuration space (see, however, Ref. for a discussion 
of alternative approaches to the constraint). These are 
important issues which will certainly deserve further in- 
vestigation. 
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Appendix A: STOCHASTIC NOISE AND LANGEVIN EQUATION 



The effective equation of motion for the in-in expectation value of the super-horizon field (p< is obtained via a 
path-integral over the sub- horizon field ip^: 



where A is the integration kernel of the action for a free scalar field, given by 



A(a;) = -a^ + 3Hdt T+m' 



(A.2) 

the effective action reads 

r[(^±] = S[ip+]-S[ip-] + i- ld^xd''x'ip'^{x)Re[U{x,x')]ip'^ix')~2 I d''xd''x'e{t-t')ip'^{x)Im[U{x,x')]ip^ {x'), (A.3) 



After some manipulations Q , setting (j)^ = Lp^ , ip^ = </i+ 



a 

and 1^*^ = (■(/)+ 



n{x,x') = dkafPt(l)kix)4'Pt'(t)l{x') , Pt = W + 3HW + 2Wdt 



where 



and the normal modes of the field arc given by (|II.5p . 

Thanks to the useful relations a^Pt — dtia^Wtik)) + 2a^Wt{k)dt and 



(A.4) 



I 2 aH (^^ 



0k = qi,{kri)4>v 



\aH) / 

after integrating by parts, the imaginary term of the effective action (|A.3p reads 



- ld%dyatai,Re /dk[^^(x)g,(A:?/) - ip'^ix)]Wtik)^k{x)Wt'ik)^lix')[ip'^{x')q*M) - ^^(x')]. 



(A.5) 



(A.6) 



Assuming Wt{k) = W{kri), then Wt = {k/at)W'{h]), where the prime denotes differentiation with respect to the 
argument of W, and with this substitution we obtain 



;.2„,k.(x-x') )H(l)*(_fc ) 

Wtik)^k{x)Wt'{k)<!>l{x') = —-^—W'{kr])W'ikr]')- 



(A.7) 



we now set 'ip^{x) = (| — v)'P<{x) + 4><i{x)/H^ and from the explicit form (|A.5p of q^, (after evaluating the integral 
over angles) we can write the integral over dk in a matrix form as 



H 



fdkk 



^sin kr W {k'q)W' [k-q') 
kr [atat')^/'^ 



Re[M^(fc?7,A:77')] 



(A., 



with the matrix M^'^ {krj, kif) given by ljII.13() . This matrix is Hermitian under the simultaneous exchange of discrete 
and continuous indices i, t — > j, t' , and Re[Mj,] is therefore symmetric. Eq. ljA.6p can thus be rewritten in the bilinear 
and symmetric form 



- ld\dV4al[ij^{x),ip^{x)]A{x,x') 



'^^{x'i 
if^ix') 



(A.9) 



where A{x,x') is the correlation matrix (|II.12|I : introducing two classical random fields and ^2 with the statistical 
weight pi. 11(1 , for the imaginary part of the effective action ((A.3p we finally get 

^-lfd-'xd\'<t>^{x)Rclnix.x')]<p^ix') ^ jj)^^J)^^p\^^^^^^y ! d^xalH[,p^{x)^^{x)+v'^(x)^^(^x)\ ^ (-^ -^q^ 

from which we immediately obtain l|II.9|l . 
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Appendix B: DISSIPATION 



111 the Langevin theory of Brownian motion, the stochastic force arising as a consequence of the colhsions with the 
particles of the thermal bath is usually split in two parts, one rapidly varying and stochastically distributed and another 
one proportional to the particle velocity, which plays the role of a friction term. The latter contribution describes 
how fast the system reaches thermal equilibrium from an out-of-equilibrium configuration, making the mean velocity 
vanish exponentially and driving the mean kinetic energy to the equipartition value. The proportionality coefficient of 
this friction term (and therefore the characteristic relaxation time needed for the system to reach equilibrium and for 
the friction to become negligible) can be derived integrating over time the correlation function of the rapidly varying 
part of the stochastic force. This relation between the macroscopic out-of-equilibrium behaviour and the microscopic 
equilibrium distribution is known as the fluctuation-dissipation theorem. 

We then expect in our situation too a dissipation term of the type aHip to show up in the equation of motion, 
with the coefficient a related in some way to the correlation functions of the noise ^. We also expect this effect to be 
negligible after a certain time. 

In Sec. ^ we disregarded the real term that appears in the action with the path-integral over the sub-horizon 
degrees of freedom (the Im[n] term in (|A.3|l '). Had we kept this term, after integrating by parts with respect to t' we 
would have obtained in the equation of motion the two extra contributions 

dya^,i}{t-t')C{x,x')ip{x') + [ d-^x'alSm^{x,x')ip{x'), (B.l) 



corresponding to non-local effects of dissipation and mass renormalization respectively, with the integration kernels 
given by 

Cix,x') ^2lTaJ dkPtMx) WrKix'), (B.2) 
5m^{x,x') = ~AS{t-t')lm [ dkWtMxWt'Ki^') ~ 2i}{t ^ t')lm ( dk Pt(f>k{x) Wt'Kix'), (B.3) 



where in Sm^ a second term proportional to S{t — t')4>]f^{x)4>'^{x') was dropped out since it gives a purely real result. 

We will concentrate on the first term of (|B.ip . showing that it is negligible compared to the usual friction term 
2>Hip\^ related to the Hubble expansion. Since it has no space dependence but e''' (^~^ going to Fourier space we 
have 



, 3 ij5 /!c4^2^'2 



[{I + ikrj) (^1 + - 2k'rj'^ (1 - tkif) 



Mt')- (B.4) 



The field (p basically contains super-horizon modes, which are slowly varying with respect to the characteristic correla- 
tion time of the sub-horizon fluctuations. We can then assume that ipk{t') is constant where the rest of the integrand 
is significantly different from zero, and take it out from the integration. We thus obtain a dissipation term ak_Hipi^(t) 
for the k mode, with a time dependent friction coefficient ak = a{h], a). Changing variables, it becomes 



2^2 



a{krj,a) — — 3— Im 



(B.5) 



this fluctuation-induced dissipation effect can then be neglected when \a{kri,a)\ ^ 3. For small values of k\ri\ 
(sufSciently after the horizon crossing of each mode), this inequality reduces to 



k\' I 



3 3-5 3-5-7 



-C 3, (B.6) 



and is always satisfied for k\'ri\ < a < 1. 

The complete behaviour of the friction cocfficcient a is shown in Figure where we see that for each mode the 
dissipative behaviour reaches its maximum for k\ri\ ^ cr, when the mode crosses the effective horizon (aH)^^; the 
spatial non-locality of the dissipation term in IjB.ip is a consequence of this effectiveness of the friction coefficient only 
about the cutoff scale for each mode 0. We also note that the sign is always negative: this "anti-dissipation" means 
for (f a global energy gain, instead of a loss, and is due to the continuous income of modes through the cutoff into 
the super- horizon field. Anyway, this effect can be neglected on a first approach, since a/3 is always much less than 
unity. 
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Figure 5: Friction coefficient a{kri, o) governing the fluctuation-induced dissipation effect for the k mode, for different values 
of (7 (from left to right, o — .125, .25, .375, .5). The sign of the coefficient is always negative (meaning energy gain), and all 
functions are peaked roughly around the cutoff scale of the mode, for k\ri\ ~ a 



Appendix C: NOISE CORRELATION MATRIX 



The noise correlation matrix Aij{x,x') can be calculated in the Gaussian window case 1)111. 6p . Since the relation 



sin(fcr)Re[M^] = -Ini[M^(e*'^'' - e^'"'')] 



(C.l) 



holds, Eq. I|II.12|I then becomes (setting S = if — r/ and R = \/ if + rj'^) 



H-6 /-oo ait^fc^ 
A(x,a;') = ^(W)'lm/ dk . 



{kri~i){kr/+i) / l+ikfj 
rpYk^ ' krj 



i{5+r)k i{5—r)k 



d 



^Im /"dfce-4^'^-'(e'(^+'-)'^- - e^^'-^^"^) 
Jo 



d 



e aR-i 



dye^ 



2 i-(5~r)cr/^/2R 

dye^ 



(C.2) 



where we have defined the derivation matrix 



d 



H 



7777' \ 00 
85 / 85^ 



1 

V 



l + i] 



8\ 8'^ 



-rji] 



86 J 86^ 
,8^ 
85^ 



(C.3) 



For r — > we have 



A(t,t')r=0 



V2H%T]f]'f I q] 8^ 



ij6 (r]rj')^ r 



e ; dye" 



R 



47r2 i?4 



(-l)'=(fc+l) (6(7 



k=0 



(2/c-l)!! \R 



2k 



( 2w_f g^(,)^^_^,^(-l)'=(fc + l)(fc + 2) fiv'-vYS'' 



(2fc- 1)! 



(C.4) 



{A^^\t — t') is given in l|III.8p 'l. and since in dc Sitter space r/ 



-Ht 



/H wc get (piL7|) . 
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Appendix D: CONDITIONAL PROBABILITY 

As a shorthand notation to skip writing ah the integration symbols, we define the inner product of two functions 
f{x) and g{x) as 

r+oo 

fg^ / d^xfix)g{x); (D.l) 



in this notation the equations of motion ljV.3p and (|V.4p become 

with (the function <t'< t") is defined to be 1 if the inequaUty is true, elsewhere) 

Ji(x') =i?(t60 < < i)^(x-x') , J2(a;') =^9(W <t' <t6o)(5(x-x') (D.3) 

Bi(.t', x") = 5(x - x')(5(x - x") / dM(i60 < t' < i)d{Un < t" < i) (DA) 

B2(t', t") = (5(x - x')5(x - x") J di§{Un <t' < i)d{Un < t" < i) (D.5) 

Since we want a probability distribution for connected correlation functions, as discussed in (jIV.4p we subtract the 
mean of the perturbative solution, which in this notation reads 

(M^]W)=-^Tr[ABi] , (J^[C](t6o)) = -^T>[AB2]. (D.6) 
If again we set Sfi = Sip — Sip^o + ■^SipQQ{t — tgo) and 6(p2 = (^Veo, we get for the joint probability (|V.l 

P{Sip,Sipeo) 



6H' 



where the new normalization A/"' reads 



1 + ^«.AB. 



-1/2 



To first order in the normalization M' and the mean value of the field cancel out, and we have; 



normalizing the previous result for the joint probability with the probability distribution (|V.2p for SfQ^, that becomes 

P{6^m) = j g-iSv2a2^-iJ^AJ2al+t^J^AB2AJ2al^ (D.IO) 

and uniforming the notation with the one adopted in Section (with jjAJj ~ Cij and A'BjAJk = Dijk), we 
obtain the conditional probability distribution l|V.5|) . 

Appendix E: COEFFICIENTS Ajfc FOR THE GAUSSIAN WINDOW 

We perform here the explicit calculation (in the small a limit) of the -Dijfc coefficients appearing in the conditional 
probability of Section expanding in the limit t — > igo and tgo ^ tin 3> . 



18 



EX — r ^ 



In 



(7^2 + fi^)2j]l„ 



1 



12 (27r)' 



-{t-teof 



(E.l) 



-D112 — -D211 



1 



" dfi {rf + r;2„ ) [f,^ + ,4 ) 2fjl^ [f,^ + ryf,, ) 



1 

12 (27r)' 



(t - teo)' (E.2) 



-D122 — ^^221 



1 



2W 4iJ ' -^^^ o In- 



77 



(77^0 + 772)27?^ 



if4 1 

T (27r)' 



-(t — t6o)(i60 ~ ^in) In 2 



(E.3) 



D121 



27ry Jb y 



In 



(l + a6)(7/ + 6) 
{y + ab)il + b) 



i?3 



4 (27r)4 



In 2 



- 1 



(E.4) 



In 



2(1 + 



(l + y)(l + 6) 



12 U' 



(E.5) 



/ 1 



^222 = — — / — 



4 „1 



8 V27r; 7fc y 



In 



(l+y)(l + 6) 



2(2/ + fo) 



(27r)^ 



H^{teo~k„f In 2 
2 



^i/^(t60 - U.? + f In^ 2 - " 



(E.6) 
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